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ABSTRACT 

The reflection nebula NGC 2023 was observed by a rocket-borne long-slit imaging spectro- 
graph in the 900 - 1400 A bandpass on 2000 February 11. A spectrum of the star, as well as that 
of the nebular scattered light, was recorded. Through the use of a Monte Carlo modeling process, 
the scattering properties of the dust were derived. The albedo is low, 0.2 - 0.4, and decreasing 
toward shorter wavelengths, while the phase function asymmetry parameter is consistent with 
highly forward-scattering grains, g ~ 0.85. The decrease in albedo, while the optical depth in- 
creases to shorter wavelengths, implies that the far-UV rise in the extinction curve is due to an 
increase in absorption efficiency 
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The far-ultraviolet absorption and scattering 
properties of interstellar dust are not well deter- 
mined, yet they are responsible for the regulation 
of the penetration of far-UV light into interstellar 
clouds. This energetic radiation can ionize atomic 
species and photodissociate molecules, and thus 
plays an important role in the ionization, chemical, 
and thermal balances of the interstellar medium. 
The extinction curve along lines of sight through 
the ISM is a measure of the combined effects of ab- 
sorption and scattering by dust, and several stud- 
ies indicate a relationship between the morphology 
of the far-UV extinction curve and diagnostics of 
the conditions within interstellar clouds (Joseph, 
Snow, & Seab 1989; Jenniskens, Ehrenfreund, & 
Desert 1992; Burgh et al. 2000). 

However, it is important to note that the ac- 
tual attenuation of radiation in a cloud depends 
strongly on the far-UV scattering properties of the 
grains responsible for the extinction, even when 
the extinction rises sharply in the far-UV (Flan- 
nery, Roberge, & Rybicki 1980). This is especially 
true for forward scattering grains. These will re- 
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move photons from line of sight observations to- 
ward stars, producing extinction, but may not ap- 
preciably affect the radiation field within clouds, 
as the photons will be only slightly deflected from 
their original direction. Thus, it is necessary to 
understand the breakdown of the extinction curve 
into absorption and scattering efficiencies sepa- 
rately, as well as the detailed properties of the dust 
itself, such as the albedo and the degree of forward 
scattering. 

Reflection nebulae provide useful sites for 
studying the scattering properties of interstellar 
dust grains. In many cases the illumination is 
due to a single star and the geometrical distribu- 
tion of the dust can be constrained, or at least 
closely approximated. As the illuminating source 
is contained very close to or within the nebula, its 
spectrum can often be measured simultaneously 
with the nebular scattered light by the same in- 
strument. This allows for the direct comparison 
of the source to the scattered light independent 
of the instrument calibration. The dust scattering 
properties may then be derived through a compar- 
ison of the observations to an appropriate model. 

A few bright nebulae have been studied in 
the ultraviolet, including the Merope nebula 
(Andriesse, Piersma, & Witt 1977; Witt 1985), 
NGC 7023 (Witt et al. 1982, 1992, 1993; Murthy 
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et al. 1993), the Scorpius OB association (Gordon 
ct al. 1994), and IC 435 (Calzctti et al. 1995). This 
work reports on the observations of the reflection 
nebula NGC 2023 with a far-UV (900-1400 A), 
rocket-borne long-slit imaging spectrographic tele- 
scope. 

NGC 2023 is a bright reflection nebula located 
- 23' SE of ( Ori. The illuminating star HD 37903 
is of B 1.5V spectral type and is embedded in the 
LI 630 molecular cloud at a distance of 450-500 pc 
(de Boer 1983). It is one of the brightest reflec- 
tion nebulae in the sky and well-studied at many 
wavelengths including [C II] 158 fim (Howe et al. 
1991; Jaffe et al. 1994), carbon recombination lines 
(Knapp et al. 1975; Wyrowski et al. 2000), far- 
infrared dust emission (Harvey, Thronson, & Gat- 
ley 1980; Mookerjea et al. 2000), and infrared H 2 
emission (Gatley ct al. 1987; Field et al. 1998; Mc- 
Cartney et al. 1999; Martini, Sellgren, & DePoy 
1999). 

A photometric study of the nebular surface 
brightness was performed by Witt, Schild, & 
Kraiman (1984, WSK hereafter) at visible wave- 
lengths. In their work, a radiative transfer model 
that incorporated multiple scattering was used to 
derive the dust scattering properties from photo- 
metric observations of the nebular surface bright- 
ness profiles in the uvby and BVRI bandpasses 
ranging from 3500 and 10000 A. In this work, 
this method is extended to far-ultraviolet wave- 
lengths; because the optical depths increase to 
shorter wavelengths, it is very important to use 
a similar detailed radiative transfer model that 
accounts for multiple scattering. 

2. Instrument and Observations 
2.1. Instrument 

The telescope flown was a modified version of 
the Faint Object Telescope (FOT) developed by 
the Johns Hopkins Sounding Rocket Group (Har- 
tig et al. 1980). The FOT comprises a 40 cm 
diameter, f/15.5 Dall-Kirkham telescope with a 
long-slit imaging spectrograph. To maximize the 
efficiency in the far-UV bandpass, the optics were 
overcoated with ion beam deposited SiC (Keski- 
Kuha et al. 1988). At the focus of the telescope is 
a mirror, which acts as the field stop for a slitjaw 
viewing camera and has cut into it a long, narrow 
(200"x 12") slit, defining the entrance aperture to 



the spectrograph. 

The spectrograph employs a 400 mm diame- 
ter Rowland Circle mounting with a = 8? 6 and 
= 0° (McCandliss et al. 1994). It contains 
a holographically-ruled, astigmatism-corrected 
diffraction grating used in —1 order. The spectro- 
graph is sealed with a motorized gatevalve, which 
operates during flight. The detector is a Z-stack 
of microchannel plates with a KBr photocath- 
ode and a delay-line anode for readout (Siegmund 
ct al. 1993). Additionally, an onboard calibra- 
tion lamp (McCandliss et al. 2000) was flown to 
provide an in-flight wavelength calibration. 

2.2. Observations 

The sounding rocket experiment was launched 
aboard a Black Brant IX sounding rocket from 
White Sands Missile Range, New Mexico (106?3 
W, 32?4 N) on 2000 February 11, at 20:27 
MST. The entire flight time was approximately 
15 minutes, with 400 seconds (from T+100 to 
T+500) spent with the detector high voltage on. 
NGC 2023 was acquired at T+175 and pointing 
stability was established through the use of a wide 
field-of-view (±1°) startracker, which maintained 
"lock" on £ Ori during the observations. 

The central star of NGC 2023, HD 37903, was 
observed for 125 seconds. However, the pointing 
was stable for only 98.5 seconds, with the star well 
situated in the slit and a steady count rate. Only 
these data are considered here and are called the 
"on-star" pointing. A far-UV spectrum of the star 
was recorded, as well as the nebular scattered light 
to a distance of ±60" along the slit, which was ori- 
ented in the north-south direction. The remainder 
of the flight was divided into three offset pointings 
to explore the scattered light distribution and to 
search for fluorescent emission from molecular hy- 
drogen (H 2 ). 

3. Data 

The data from the rocket experiment were 
telemetered to the ground in the form of a bit 
stream. For every detected photon, as well as cal- 
ibration lamp photons, dark counts, and possible 
detector hot spots, X, Y, and T information was 
encoded in the stream. X and Y correspond to 
the position on the detector of the event and T 
the time of the event during flight. These data 
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may be binned spatially and/or temporally. 

A spectral image was created, shown in Figure 

1, by compressing the entire flight's data into an 
X-Y image. The image has been corrected for sev- 
eral detector effects, as well as for spacecraft jitter 
during the "on-star" portion. Present in the image 
are the continuum flux from HD 37903 (the hori- 
zontal stripe) and the nebular scatter surrounding- 
it. The strong vertical line at 1216 A is the H I 
Ly-a airglow emission. Airglow lines at 1026 A 
and 1304 A are also present. At the top of the 
image is a spectrum of atomic, ionic, and molec- 
ular nitrogen provided by the onboard calibration 
lamp. As the light from the lamp did not enter the 
spectrograph through the entrance slit, there is a 
slight zero-point offset in the dispersion direction. 

The calibrated stellar flux, computed from the 
central 10 pixels of the detector, is shown in Figure 

2. The flux is consistent with measurements made 
of the star by the Hopkins Ultraviolet Telescope 
(HUT) (Buss ct al. 1994). Clearly seen in the spec- 
trum are several interstellar and stellar absorption 
features, including those of H I, C II, C III, O I, 
Si II, and the H 2 Lyman and Werner bands. Also, 
the strong effect of dust extinction is apparent, 
since the intrinsic spectral energy distribution of 
a B1.5V star rises to shorter wavelengths. The fea- 
ture marked "grid" is an instrumental artifact due 
to the vignetting of light by the QE enhancement 
grid on the detector. 

Due to the two-dimensional nature of the long- 
slit imaging spectrograph, a spatial profile of the 
nebular scattered light can also be extracted from 
the data. The data were compressed in the spec- 
tral dimension, excluding the region around Ly- 
a, binned up spatially along the slit as a surface 
brightness, assuming that the scattered light fills 
the slit, and referenced to the flux of the star. This 
produced a slit profile of the ratio of the nebular 
surface brightness to stellar flux (S/F*) (see Fig- 
ure 3), which is independent of the absolute detec- 
tion efficiency of the spectrograph. Ovcrplotted 
is the instrumental line spread function (LSF) as 
determined from laboratory measurements using a 
windowless vacuum ultraviolet collimator (Burgh 
et al. 2001). The inner region is dominated by 
light from the star, but outside of ±10" there is 
negligible contribution from the LSF and the neb- 
ular flux can clearly be observed, falling off expo- 
nentially with distance from the central star. To 



the south of the central star, a dust lane appears 
to intersect the nebulosity, causing a depression in 
the flux centered at ~ 50" south. 

4. Extinction Curve and Ry 

The derivation of dust scattering properties in 
reflection nebulae is strongly model dependent. 
Relevant parameters include the physical depth of 
embedment of the illuminating star in the dust 
cloud and the wavelength dependence of the dust 
opacity. A depth of 1 pc was assumed by WSK. 
This is consistent with the size of the nebula as 
determined from the emission measure of carbon 
recombination lines (Knapp et al. 1975). At the 
distance of the L1630 cloud, this corresponds to 
<~ 450". The wavelength dependence of the dust 
optical depth is derived from the extinction curve; 
however, it is very important to determine what 
fraction of the extinction along the line of sight to 
the illuminating star is associated with the nebula 
itself and what fraction is due to foreground ma- 
terial. The optical extinction through the L1630 
cloud is estimated at Ay ~ 40 mag (Harvey, 
Thronson, & Gatley 1980); because the visible ex- 
tinction toward HD 37903 is Ay ~ 1.4 mag, it is 
clear that the star is located on the near side of this 
dense cloud. From far-infrared (FIR) measure- 
ments, Harvey, Thronson, & Gatley (1980) sug- 
gest that the extinction from FIR-cmitting dust 
should account for a significant fraction of the to- 
tal extinction. Since this dust emission is due to 
the processing of UV photons, it is likely that most 
of the extinction along the line of sight is associ- 
ated with NGC 2023. 

To determine the extinction curve, the "pair 
method" was used, in which the spectrum of an 
unreddened star of the same spectral type is com- 
pared to that of HD 37903. Following de Boer 
(1983), HD 37744 was used as the comparison star 
for HD 37903. This star, only <~ 37' away from 
HD 37903, lies to the west of the L1630 cloud and 
exhibits only a mild amount of reddening (Eb-v 
= 0.05) (Hardie, Heiser, & Tolbert 1964; Craw- 
ford ct al. 1971). Its spectroscopic parallax sug- 
gests that it lies only ~ 50 pc more distant than 
HD 37903 (de Boer 1983), and as such samples 
essentially the same foreground extinction. The 
difference in reddening is thus due to the dust in 
the LI 630 cloud that lies between the front surface 
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of the cloud and the embedded star, and the ex- 
tinction curve derived from comparing these two 
stars should be indicative of the extinction prop- 
erties intrinsic to the nebular dust, canceling out 
the effects of the foreground dust along the Orion 
line of sight. 

Low-dispersion data of both HD 37903 and 
HD 37744 were obtained from the International 
Ultraviolet Explorer (IUE) archive. For each star, 
the short- and the long-wavelength data were 
joined to a common wavelength grid ranging from 
1150 to 3350 A. The extinction curve, shown in 
Figure 4, was derived from the ratio of the stel- 
lar fluxes, normalized to the difference in E B _y, 
and fitted with the parameterization of Fitzpatrick 
& Massa (1990). The extinction curve, k(x) = 
E\_ v I E B _ V with x = A -1 , is converted to opti- 
cal depth through the use of the ratio of total-to- 
selective extinction, Ry [= A v /E B _ V ], following 
r = (k(x) + R v )E B -v /1.086. 

The value for Ry was reassessed for this work. 
The extrapolation of the extinction curve to infi- 
nite wavelength, or zero frequency, at which ab- 
solute extinction (A\) would be zero, provides for 
the calculation of Ry. The value of Ry thus de- 
termined depends strongly on the functional form 
of the extinction toward zero frequency. It is use- 
ful to have extinction data at low frequency and 
near-infrared (NIR) photometric observations are 
typically used. 

It has been found that the shape of the NIR 
extinction curve is very uniform along almost all 
lines of sight (Rieke & Lebofsky 1985; Whittet 
1988; Cardelli, Clayton, & Mathis 1989; Martin 
and Whittet 1990), taking the form of a power law 
(A\ cx A~") between ~ 1 ^m and 5 /im. At longer 
wavelengths, emission from dust can contaminate 
photometric observations and complicate the de- 
termination of extinction. Martin and Whittet 
(1990) expand this power law in terms of the ex- 
tinction curve and Ry in the following form: 

= eA"<* - i^- (1) 

For HD 37903, the NIR photometric data from 
the literature, listed in Table 1, were averaged to- 
gether and the colors A — V determined. The pho- 
tometry for V and B — V are taken from Hardie, 
Heiser, & Tolbcrt (1964) (V=7.8A, B-V=0A0), 
Lee (1968) (V=7.83, 5-^=0.11), and Racine 



(1968) (V=7.82, B-V=0.09). These were av- 
eraged together to produce a V magnitude and 
B — V color. By comparing the NIR colors to 
the intrinsic colors for a B1.5V star (Koornneef 
1983), the color excesses E\_y were computed. 
These were then normalized by E B -y and a least 
squares fit was performed with the Martin and 
Whittet (1990) power law described above (see 
Figure 5). This method produced a determina- 
tion of a = 1.8 ± 0.9 and R v = 3.83 ± 0.40. The 
large error in these values is mainly due to the 
uncertainty in the L and M band photometry. 

These results compare well with previous deter- 
minations of a and Ry. Whittet (1988), from an 
analysis of an extensive data set obtained from 
the available literature at the time, found an 
a = 1.70 ± .08. Cardelli, Clayton, & Mathis 
(1989) used the Rieke & Lebofsky (1985) curve, 
which has an a ~ 1.6 and Landini et al. (1984), 
through observations of the hydrogen recombina- 
tion lines in H II regions, found values for a and 
Ry of 1.85 ± 0.05 and 4.11 respectively. Ry can 
also be computed from color excesses following 
Ry = lAEy_j(/E B _y, from the standard van 
dcr Hulst theoretical curve, and gives a value of 
Ry = 3.7. These are comparable to the value of 
Ry based on the above calculation within the error 
introduced by the uncertainty in the NIR photo- 
metric measurements. 

5. Monte Carlo Model Fitting 

An analytic solution for the surface brightness 
of a reflection nebula can be computed in the case 
of single scattering. However, this would only be 
accurate for optically thin nebulae. Roark et al. 
(1974) have shown, through the use of Monte 
Carlo methods, that multiple scattering can have 
a significant effect on the observed surface bright- 
ness distribution. The intrinsic random nature 
and statistical approach to the scattering process 
has allowed for radiative transfer models of dusty 
nebulae to be constructed, including some with 
arbitrary dust density and source distributions 
(Witt 1977; Yusef-Zadeh, Morris, & White 1984; 
Wood and Reynolds 1999; Gordon et al. 2001). 

To determine the dust scattering properties in 
NGC 2023, a Monte Carlo dust radiative transfer 
model was developed following the algorithm de- 
scribed in full detail by Gordon et al. (2001). To 
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simplify the modeling process and reduce the com- 
putational demands, the nebula was modeled as a 
single star embedded in a spherical dust cloud of 
constant density and finite extent. 

Using this Monte Carlo process, a model nebula 
can be created for a given radial optical depth, an- 
gular size of the nebula, and the dust's albedo, a, 
and phase function asymmetry parameter, g, as- 
suming the phase function of Henyey & Greenstein 
(1941). The first two are set and the dust prop- 
erties a and g are varied until agreement with the 
observations are maximized. The following section 
discusses the results of fitting models to the data 
and the resultant derivation of the dust scattering 
properties. 

5.1. Radial Profile Fitting 

The data were binned spatially to produce a slit 
profile, compressing the entire wavelength range 
into a one-dimensional profile, excluding a 10 A 
band (one slit width) around the Ly-a geocoronal 
airglow line at 1216 A. The error for each bin was 
computed using Poisson statistics, i.e. Uhin — VN~, 
with TV the number of counts in the spatial bin. 
The inner ±5'.'5 of the profile is dominated by the 
star and we take the total counts in this region to 
correspond to the stellar flux (F»). 

The flat- field corrected ratio of the nebular sur- 
face brightness to stellar flux, S/F*, for the "on- 
star" exposure is shown with logarithmic scaling in 
Figure 3. Positive offsets are to the north. Outside 
of the inner region, where the instrumental LSF 
causes the stellar flux to dominate, the nebular 
brightness falls off exponentially by about an or- 
der of magnitude. Only the region inside of ±60" 
contained data with a count rate above the dark 
rate. To the south, a dust lane cuts in front of the 
nebulosity at around —50". 

A grid of models for various values of a and 
g corresponding to the far-UV optical depth was 
generated. The optical depth used was a stellar 
flux intensity weighted average of the t(A) derived 
from the extinction curve in Section 4. For each 
model in the grid, the \ 2 statistic was computed 
by comparing the model to the data following 

X 2 =Y^ y ( Xi ' , T,a,9 " > ^j 2 (2) 

In this definition, yi and cr^ are the data and their 



respective errors for the ith spatial bin. The term 
y(xi;T,a,g) is the model prediction for S/F* as a 
function of radius for the given optical depth and 
the input values for a and g. The sum is over the 
region r = 10 — 60" on both the north and the 
south sides of the central star. The instrumen- 
tal LSF was added in to the model fits; however, 
its presence had no effect on the fitted parameters 
as the contribution from the star due to the LSF 
is negligible outside of 10". The best fit param- 
eters were determined from the model that gave 
the minimum value of x 2 . The best fit parameters 
for the full experiment bandpass are listed at the 
top of Table 2 and the model is shown overplotted 
on the data in Figure 6. Additionally, as the pro- 
file showed some asymmetry, the north and south 
sides of the star were fitted separately. These pro- 
duced different fit results; however, that may be 
more likely due to a variation in the dust distribu- 
tion than actual differences in the scattering char- 
acteristics. 

To obtain the wavelength dependence of the 
dust scattering properties, the data were divided 
into several wavelength regions. Six wavelength 
regions were chosen across the full bandpass, and 
the above procedure was used to fit these data, 
with a few exceptions. The spatial range was lim- 
ited to less than 60". This is due to the fact that 
as fewer data were used the contrast between the 
nebular scattered light and the background was 
worse and the outer limit of the range was de- 
creased to avoid fitting the background. The max- 
imum fitting radius was chosen for each bandpass 
so that the signal-to-noise ratio (S/N) of the spa- 
tial bins was always greater than one. The widths 
of the bandpasses were chosen so that this ra- 
dius was greater than 30". The wavelength ranges 
and maximum radii are listed in Table 2. For the 
shorter wavelength bandpasses, in which there are 
considerably less data due to the linc-of-sight ex- 
tinction, a coarser spatial binning was used. This 
provides more counts per spatial bin, which results 
in a higher S/N per bin, and thus a larger radius 
before the S/N dropped to one. If the maximum 
radius were too close to the star, the small num- 
ber of data points would limit the constraint on g. 
The best fit dust parameters are listed in Table 2. 

Confidence intervals were based on the x 2 con- 
tours. Figure 7 shows the \ 2 contours for the 
1240-1295 A bandpass. The x 2 fitting produces 
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asymmetric error bars, with most of the confidence 
tending toward lower g and higher a. Thus, it is 
apparent that there is a strong coupling between 
a and g. Any random or systematic effects that 
would change the value of g could have a corre- 
sponding change in a. 

The derived dust parameters from Table 2 are 
plotted in Figure 8. For comparison, the predicted 
values from Weingartner & Draine (2001) for val- 
ues of Ry of 3.1, 4.0, and 5.5 are overplotted. We 
find g increases mildly across the bandpass, while 
a falls toward shorter wavelengths. 

5.2. Spectral Fitting 

To achieve a better S/N, the data can be 
binned up spatially and the ratio of the nebular 
surface brightness to stellar flux shown as a func- 
tion of wavelength rather than slit position. How- 
ever, this removes the ability to fit g in an inde- 
pendent manner, so for this procedure a constant 
value of g was assumed, as supported by the re- 
sults of the previous section. Then, as a function 
of wavelength, the S/ F* can be used to derive the 
albedo directly through the model fitting. 

The data are binned in the spectral dimension, 
compressing the spatial range from 10 - 60" into 
one bin for each wavelength range. Because the 
nebular scattered light is a diffuse source and is 
expected to fill the slit, the data were binned by 
one slit width, or 10 A across the full wavelength 
region, excluding the Ly-a region, to 1030 A. 
At shorter wavelengths, the detected count rate 
for the nebular scattered light is consistent with 
the background. The stellar spectrum was also 
smoothed by one slit width before the ratio was 
determined. The derived S/F* is shown in Fig- 
ure 9. It is virtually constant at all wavelengths, 
despite a rising optical depth toward shorter wave- 
lengths. 

The data arc compared to a grid of models with 
different albedos for an assumed value of g = 0.85 
across the bandpass. Since the g from the ra- 
dial profile fitting was fairly constant, the average 
value is used in the fitting here. The best fit albe- 
dos are shown in Figure 10. There is a decrease in 
albedo as the wavelength decreases. This is con- 
sistent with the results derived from the profile 
fitting above (see Figure 8). Overplotted are the 
predictions from Weingartner & Draine (2001) for 



a few values of Ry. The derived albedos agree 
very well with the predictions. 

The relative efficiencies for both scattering and 
absorption, Q S catt and <5 a b s , can be derived from 
the albedo, since a — Q SC att/Qcxt- Therefore, 
Qscatt = aQoxt and <3 abs = (1 - a)Q cyLt . The ex- 
tinction efficiency, Q oxt (= Q a b s + Qscatt), is pro- 
portional to the extinction curve, and thus to r. 
The derived relative efficiencies are shown in Fig- 
ure 11, along with the extinction curve, in units 
of Ay/Es-v- It is apparent that the rise in the 
far-UV extinction curve of NGC 2023 is due to 
an increase in absorption, whereas the scattering 
efficiency is nearly constant. 

5.3. Systematics 

As mentioned earlier, the derivation of the dust 
parameters a and g is sensitive to the modeling 
process. There are several sources of systematic 
error that can affect the best fit parameters. This 
section discusses some of these sources and de- 
scribes the effects they have on the derived dust 
properties. 

In the data preparation process, it was assumed 
that all of the detected photons within ±5'/5 of 
the peak of the radial profile are stellar in origin. 
However, a contribution from scattered light at 
small angles, which would be included with the 
stellar light, would result in an underestimation of 
S/F* and hence the derived albedo. This effect 
would become increasingly important as all the 
parameters, r, a, and g, increase. As the optical 
depth increases, the contrast between the escaping 
stellar light, F*, and the nebular scattered light 
would decrease due to the increased probability of 
interaction. An increase in albedo would make the 
nebular scattered light brighter, and higher values 
of g increase the chance of scattering into the line 
of sight at small angles from the star. 

Assuming that the constant density of the 
model can be continued to the position of the 
star, the Monte Carlo model was used to estimate 
the number of nebular photons near the center 
that would be mixed in with the stellar light. For 
r = 2 - 3, a = 0.2 - 0.4, and g = 0.8 - 0.9, 
the fraction of nebular contribution to the stel- 
lar flux ranges from 5 - 20%. However, based on 
C II emission, Knapp et al. (1975) conclude that 
there may be a dust-free H II region surrounding 
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HD 37903 with a radius as large as ~ 25", and 
Harvey, Thronson, & Gatley (1980) suggest that 
the inner 0.04 pc (18") is clear of dust. A "hole" in 
the dust distribution would reduce the contribu- 
tion from scattered light in the vicinity of the star, 
and the modeled contamination is thus an upper 
limit. Decreasing the stellar flux in the model 
fits as described above resulted in an increase in 
the fitted albedo by an amount proportional to 
the flux reduction, but did not affect <?; i.e., a de- 
crease in flux by 10% corresponded to an increase 
in the derived value of a by 10%. 

A spherical geometry was assumed in the mod- 
eling process. For the L1630 cloud, a plane- 
parallel slab model may be more appropriate. 
However, if the depth of embedment, and thus 
the radius of the sphere in the model, is large rel- 
ative to the extent being observed as projected 
onto the sky, the difference between the sphere 
and slab models is negligible. We used the value 
from WSK of 1 pc for the radius, and therefore the 
model fits were performed on the inner 0.5—1 '0 of 
an assumed ~ 7' nebula. If the depth of embed- 
ment, and hence the radius of the spherical model, 
is of lower value, then the fitted g values decrease. 
Also, the correlation of a and g, as apparent from 
the x 2 confidence intervals shown in Figure 7, is 
weaker for smaller depth. 

At too small a radius relative to the observed 
extent, the assumption of spherical symmetry is no 
longer sufficient to describe the geometry. In this 
case, it is possible that the observed nebular sur- 
face brightness distribution can be explained by a 
slab geometry of shallow depth but large horizon- 
tal extent. There do exist degeneracies between 
the two geometries, such that a spherical cloud of 
given a and g can appear identical to a slab with 
the star embedded less deep and lower values of 
g but observed out to radii larger than the actual 
depth of embedment. The value of 1 pc for the 
radius produces fits to the visible wavelength sur- 
face brightness with reasonable values for a and g 
(WSK). 

The optical depth used in the modeling is de- 
pendent on the far-UV extinction curve. The pri- 
mary systematic error affecting t\jv is from the 
propagation of the L and M band photometric er- 
rors in the determination of Although these 
were large, the resultant uncertainty in Ry does 
not have a large effect on the optical depth in the 



UV and the fitted parameters are modified only 
slightly (< 10%). If less of the extinction were as- 
sociated with the nebula itself, the albedo would 
increase; however, not more than half can be fore- 
ground or the far-UV albedo would be driven to 
the maximum value of 1.0. There is also an im- 
plicit assumption that, for the lack of a far-UV 
spectrum of HD 37744, the UV extinction curve as 
derived from IUE can be extrapolated to shorter 
wavelengths than 1150 A following the parameter- 
ization of Fitzpatrick & Massa (1990). 

As already mentioned, a constant dust density 
distribution was assumed in the modeling process. 
This greatly reduced the computation time, be- 
cause the distance each photon traveled was sim- 
ply linearly proportional to the optical depth along 
all directions through the model dust cloud. WSK, 
on the other hand, assumed a weak radial falloff in 
the density distribution. Such a radial density de- 
pendence could account for some of the observed 
falloff in nebular surface brightness and reduce the 
value of g necessary to fit the slit profiles. Al- 
though there is the possibility of large scale density 
variations, the nebula appears to be fairly smooth 
in both the blue, as observed by the digitized sky 
survey, and the far-UV, as shown in the slit profile 
in Figure 3. 

It should be noted that the predicted dust prop- 
erties, such as those of Weingartner & Drainc 
(2001), to which the results have been compared, 
are based on Mie scattering theory, which solves 
for the exact electromagnetic fields in a spheri- 
cal or cylindrical dust grain. However, it is clear 
from polarization studies that dust grains in the 
interstellar medium are not spherical. The more 
probable situation is that dust grains collide and 
either stick together or shatter, producing grains 
of a more random distribution of shapes (Mathis 
& Whiffen 1989; Bazell & Dwek 1990; Fogel & Le- 
ung 1998, e.g.). Savage & Mathis (1979) point 
out that the Henyey-Greenstein phase function is 
an analytic function that approximates the results 
from Mie theory and may not represent the true 
scattering phase function. 

Witt et al. (1993) have pointed out that dust 
scattering is often accompanied by fluorescent 
emission of H2. It is therefore necessary to ac- 
count for any contribution from the H 2 far-UV 
fluorescent cascade to the dust scattered contin- 
uum. Fluorescently excited H 2 is indeed present 
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as evidenced indirectly by the detection of absorp- 
tion lines from highly excited vibrational states 
(Meyer et al. 2001). More direct evidence is of- 
fered by the IUE spectrum, SWP 38082, which 
exhibits the same characteristic pair of fluorese- 
cent H 2 emission lines at 1577 and 1608 A as seen 
in the IC 63 (SWP 33898) spectrum where H 2 flu- 
orescence has been confirmed (Witt et al. 1989). 
The scattered continuum in NGC 2023 has an 
S/F* ~ 10 5 ' 6 in the wavelength region longward 
of 1700 A, and the H 2 emission peaks at <~ 3 times 
this value. 

However, the IUE spectrum of NGC 2023 was 
acquired at a position <~ 75" east-north-east of 
HD 37903, coincident with a bright infrared emit- 
ting knot as revealed in high spatial resolution im- 
ages of the nebula. These images, which isolate 
specific vibrational quadrupolar transitions in the 
ground electronic state of H2, generally show a 
clumpy distribution (Field et al. 1998; Takami et 
al. 2000). In particular, the Takami et al. (2000) 
images shows the region to the north of HD 37903 
to be devoid of H 2 infrared emission within the ex- 
tent of the rocket experiment's slit. To the south 
there are two lanes; one is rather close, ~ 10 - 20", 
to the central star and comparable in brightness 
to the east-north-east knot, and the other (the 
southern bar) is - 70" south of HD 37903 and 
brighter, with a void of emission in between. Mar- 
tini, Sellgren, & DePoy (1999) have observed the 
infrared emission spectroscopically in the region 
to the south of the star that is spatially coinci- 
dent with our slit. They find that the brightest 
infrared H 2 emission lies in a bar between 67" and 
84" south that is <~ 3.5 times stronger than the 
region closer to the star. They also show that the 
infrared continuum emission is essentially negli- 
gible in the southern bar and increases towards 
the star, becoming comparable to the H 2 emission 
close in. 

The southern bar of the nebula was observed 
within our slit; however, the data here are consis- 
tent with the background, with a 2 a upper limit 
of S/F* — 10 5 5 . As this is the brightest region 
of H 2 infrared emission along the slit, we conclude 
that the region closer to the star, in which we de- 
tect S/F* = 10 6 - 5 — 10 7 , has a neglible contribution 
from far-UV H 2 fluorescent emission. The asym- 
metry between the profiles in the north and the 
south sides of the star mentioned previously is thus 



more likely due to asymmetries in the dust distri- 
bution than additional emission from H 2 . This 
asymmetry has only a small effect on the fitting of 
the modeled slit profile to the data. 

We note that the rocket flight included two off- 
set pointings in the north-eastern portion of the 
nebula, ^45" from the star, which do show an ex- 
cess of emission from the best-fit scattered light 
model by a factor of <~ 2 — 4. Although the spec- 
tral resolution of the sounding rocket is ideal for 
investigating dust scattered light it is insufficient 
for unambiguously detecting suspected far-UV flu- 
orescent emissions from H 2 . 

6. Discussion and Conclusions 

The results described above show that the dust 
in the reflection nebula NGC 2023 is highly for- 
ward scattering and has a low albedo that de- 
creases with decreasing wavelength. These results 
are in good agreement with the predictions of a 
and g based on the fitting of extinction curves (Li 
k Draine 2001; Weingartner & Drainc 2001). The 
rise in the far-UV extinction curve is due to an 
increase in the absorption efficiency of the dust 
grains, whereas the scattering efficiency appears 
constant. This suggests that the scattering is be- 
ing performed by large grains, while a low-albedo 
small grain population may be the cause of the 
short wavelength absorption. 

The approach taken here is very similar to 
that of previous work on other reflection nebulae. 
NGC 7023, for example, is a similar, but more ex- 
tincted, nebula to NGC 2023. It has been observed 
extensively in the ultraviolet, by IUE (Witt et al. 
1982), the Ultraviolet Imaging Telescope (UIT) 
(Witt et al. 1992) and the Hopkins Ultraviolet 
Telescope (HUT) (Murthy et al. 1993) as part of 
the Astro- 1 mission, and Voyager 2 (Witt et al. 
1993). 

The initial study of NGC 7023 with IUE sug- 
gested more isotropic scattering in the UV than 
the visible, but the measurements with UIT pro- 
vided more detailed surface brightness profiles 
and Witt et al. (1992) conclude that the dust 
in that nebula is forward scattering (g ~ 0.70). 
The UIT measurements indicated a high albedo 
(a <~ 0.6), similar to that in the visible, but the 
short wavelength limit was 1400 A. The HUT mea- 
surements extended into the windowless ultravio- 
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let bandpass covered by the sounding rocket ex- 
periment. Murthy et al. (1993), using a similar 
Monte Carlo method, also found that the albedo 
was high (~ 0.6 - 0.7) in the 1000-1400 A range, 
but decreasing to shorter wavelengths. They note 
that their derived albedo would decrease if more of 
the dust can be associated with the nebula. Witt 
ct al. (1993) show models for three different op- 
tical depth radii in their fitting of the Voyager 
data, resulting in dust albedos at A = 1000 A of 
0.45 - 0.26, with each case showing a clear decline 
(about 25%) when going from 1400 to 1000 A. 

Calzetti et al. (1995) observed the reflection 
nebula IC 435, which is very near NGC 2023 in 
Orion. Their data were taken with IUE and were 
intended to investigate the scattering properties 
near the 2200 A bump. Also employing a Monte 
Carlo model, they derived the dust albedo and 
phase function asymmetry parameter. The results 
indicate a drop in albedo at the bump, but rising 
albedo toward the short wavelength limit of their 
observations. 

It might be expected that, because this neb- 
ula appears to also be in the L1630 cloud, the 
dust characteristics would be similar to those of 
NGC 2023. However, HD 38087, the central star of 
IC 435, appears to have a lower than average far- 
UV rise, and Cardelli, Clayton, & Mathis (1989) 
find an Ry of 5.3. Aiello et al. (1988) list a value of 
R v = 4.86 for HD 38087, but it is clear that the 
extinction toward this nebula is more consistent 
with a dense region where the dust grains may be 
larger than average. Also, as this star is of spectral 
type B5V, the dust environment around the star 
may be different than that of NGC 2023 if any pro- 
cessing of the dust grains occurs due to the high 
ultraviolet flux of the hotter HD 37903. However, 
Calzetti et al. associate only ~ 20% of the ex- 
tinction along the line of sight to the nebula itself. 
This would indicate that the foreground material, 
which would therefore account for most of the ex- 
tinction curve, has the dust properties of a dense 
environment. This is inconsistent with the calcu- 
lation of the extinction curve toward HD 37903 as 
described in Section 4, which indicates little fore- 
ground extinction toward the L1630 cloud. Had 
they used a higher optical depth in their Monte 
Carlo model, their derived albedos would be of 
lower value, more in agreement with the results 
derived here. 



Measurements of dust scattering properties 
through observations of the diffuse galactic light 
(DGL) have significant scatter. Some recent ultra- 
violet measurements have been made. Observa- 
tions with the Berkeley UVX spectrometer (Hur- 
witz, Bowyer, & Martin 1991) indicate a low a, 
but also a low g in the 1400-1850 A bandpass. 
Witt, Fricdmann, & Sasscen (1997) found, with 
the Far Ultraviolet Space Telescope (FAUST) on 
board the Atlas- 1 mission, an albedo of 0.45 and 
g = 0.68 in the same bandpass. Schiminovich, 
Friedman, Martin, & Morrissey (2001) obtained 
similar results from modeling the observations 
from a wide-field ultraviolet sounding rocket ex- 
periment, finding a = 0.45±0.05 and g = 0.77±0.1 
at 1740 A. These slightly lower albedos for the 
diffuse ISM, relative to those found in the denser 
regions of NGC 7023 and IC 435, have been asso- 
ciated with a dust size distribution including more 
small grains than average. 

There is a large spatial variability to Galactic 
interstellar extinction, particularly in the ultravi- 
olet. Witt et al. (1993) point out that if these 
variations are due to changes in the contribution 
of various grain components to the total extinc- 
tion then it may be expected that nebula in differ- 
ing environments may show different far-UV dust 
albedos. 

Though agreeing with the predictions of Wcin- 
gartner & Draine (2001), the albedo of the dust 
in NGC 2023 appears to be more consistent with 
the observations of lower albedo dust in the dif- 
fuse ISM as determined from measurements of the 
DGL, even though it resides within a reflection 
nebula with an Ry indicating a denser than aver- 
age environment. Some of this discrepancy could 
be due to the assumption that the dust in the 
nebula is homogeneous. Witt & Gordon (1996) 
have performed simulations of the radiative trans- 
fer through two-phase clumpy dust of various fill- 
ing fractions. Their models suggest that radiative 
transfer techniques based on homogeneous mod- 
els can underestimate the albedo when applied 
to inhomogeneous nebulae. This effect becomes 
more significant at higher values of g and r. In 
their study of the infrared scattering properties in 
NGC 2023, Sellgren, Werner, & Dinerstein (1992) 
suggest that the dust may be clumpy in order to 
reconcile the J, H, and K band surface brightness 
measurements with the UV and visible measurc- 
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ments. 

It is critical to determine the far-UV scatter- 
ing properties of dust in order to understand the 
depth dependent effects of the penetration of far- 
UV photons into interstellar clouds. In their 
study of the photodissociation of carbon-bearing 
molecules, van Dishoeck & Black (1989, vDB here- 
after) show that varying the morphology of the 
extinction curve can cause order of magnitude 
changes in the central photorates in translucent 
clouds, in particular for CN and CO, which pho- 
todissociate at wavelengths shorter than 1100 A. 
These effects may have been observed in the stud- 
ies of Joseph, Snow, & Seab (1989) and Burgh 
et al. (2000). However, it should be noted that the 
vDB study assumed the dust scattering character- 
istics of Model 2 of Roberge, Dalgarno, & Flan- 
nery (1981), which has an a = 0.6 and g = 0.5. 
These values are not consistent with the results 
found here, nor with many of the other studies, 
nearly all of which find a higher value of g in 
the ultraviolet. Had vDB used the dust proper- 
ties indicated in this study of NGC 2023, their 
photorates may have been different. The lower 
albedos indicated here would drop the dissocia- 
tion rates due to the stronger attenuation by the 
dust; however, the more forward-scattering grains 
could possibly counteract that effect by increas- 
ing the probability that a photon might get deeper 
into the cloud before interacting with either a dust 
grain or molecule. 

Although it is very important to determine the 
ultraviolet scattering properties of dust, it is clear 
that there are several inherent difficulties in the 
modeling of scattered light in reflection nebulae. 
The most important being the geometry and dis- 
tribution of the dust relative to the illuminating 
source. Further studies of the far-infrared emission 
from dust, in conjunction with these and other far- 
ultraviolet measurements of dust scattered light, 
may help to constrain the geometry and thus limit 
the introduction of systematic errors into the de- 
termination of the dust scattering properties. 
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Table 1 

Near-infrared photometry of 
HD 37903. 



Magnitude 



Band 


A (/xm) 


Ref 1 


Ref 2 


Ref 3 


I 


0.88 


7.52 






J 


1.25 


7.48 


7.36 


7.44 


H 


1.65 




7.28 


7.35 


K 


2.20 


7.41 


7.28 


7.32 


L 


3.50 


7.00 


7.44 


7.36 


M 


4.80 






7.20 



References.— Ref 1 - Lee (1968); Ref 2 
- Whittet & van Breda (1978); Ref 3 - DePoy 
ct al. (1990) 



Table 2 

Best fit a and g from radial profile fitting. 



Bandpass 



A m in (A) A max (A) 



Max.* 
radius (") 



Best-fit dust parameters 
a g 



990 


1380 


60 


2.8 


0.26 


0.89 


990 


1380 


60 N 


2.8 


0.27 


0.87 


990 


1380 


40 S 


2.8 


0.23 


0.92 


1050 


1115 


33 


3.6 


n 9 o+0.20 
u - zo -0.05 


o-88±8:2S 


1115 


1150 


33 


3.2 


21 +0 - 25 

u - zl -0.05 


o.89l8:S 


1150 


1190 


35 


3.0 


n r, 9 +0.25 
u - zz -0.03 


90+ 04 
u.au_ 23 


1240 


1295 


40 


2.7 


S1+ 012 


0.851°;?^ 


1310 


1345 


40 


2.5 


o.35±g:ig 


0.83i°;l° 


1345 


1380 


35 


2.4 


0-39t°; 2 


o.so±ns 



t Maximum radius of the fitting region 

* Intensity weighted optical depth for the given bandpass 
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Fig. 1.— Long-slit spectral image of NGC 2023. 
The long axis of the slit was oriented with north at 
the top and south at the bottom. The spectrum 
at the top is that of the onboard calibration lamp. 
Red indicates the highest count rate regions. 
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Fig. 4. — Derived extinction curve for HD 37903 
from IUE data. The solid line is a fit following the 
parameterization of Fitzpatrick & Massa (1990). 
The dashed line is the extinction curve following 
the formula of Cardelli, Clayton, & Mathis (1989) 
for R v = 3.83. 



Fig. 2. — Flux calibrated spectrum of HD 37903. 
The error in the flux is shown as the dashed 
line. Strong stellar and interstellar absorptions 
are identified. The feature marked "grid" is an 
instrumental artifact. 
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Fig. 3. — Spatial profile of nebular scattered light. 
Filled circles with error bars are the data, the solid 
line is the instrumental line spread function (LSF). 
The profile is centered on HD 37903 with positive 
offsets to the north. 
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Fig. 5. — Infrared color excesses with extinction 
fit overplotted. Shown are the best fit parame- 
ters following the offset power law of Martin and 
Whittet (1990). 
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Fig. 6. — Slit profile of nebular surface bright- 
ness to stellar flux ratio for the "on-star" position. 
These data encompass the entire detector band- 
pass, excluding Ly-a. Overplotted is the model 
corresponding to the best fit a and g. The profile 
is centered on HD 37903 with positive offsets to 
the north. 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 L 
0.0 




0.2 



1 o (68%) 

2 a (90%) 

j i , , , i i 

0.4 0.6 
a 



0.8 



1.0 



Fig. 7. — Confidence intervals from the \ 2 cori ~ 
tours for the model fit of the nebular brightness 
radial profile for the 1240-1295 A bandpass. The 
black dot indicates the best fit parameter values, 
corresponding to the minimum \ 2 . There is a de- 
generacy between the two parameters and the con- 
fidence tends toward lower g and higher a. The 
confidence intervals for the other bandpass fits are 
essentially of similar shape. 
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Fig. 8. — Derived values of a and g for the profile 
fits listed in Table 2. Theoretical predictions for 
various values of Ry from Wcingartner & Drainc 
(2001) are overplotted. 
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Fig. 9. — Spectrally binned S/F*. The region 
around Ly-a is excluded. 
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Fig. 10. — Derived wavelength dependence of 
albedo from model fitting of the spectrally binned 
data. The prediction for Ry values of 3.1, 4.0, and 
5.5 from Weingartner & Draine (2001) are over- 
plotted. 
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Fig. 11. — Relative values of the extinction, ab- 
sorption, and scattering efficiencies derived from 
the best fit albedos. The dust extinction efficiency, 
Q ext , is proportional to the extinction curve for 
HD 37903. 
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